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Abstract 



. We present numerical evidence for the existence of new black hole solutions in d > 6 spacetime 

dimensions. They approach asymptotically the Minkowski background and have an event horizon topol- 
ogy X 5"'"''. These static solutions share the basic properties of the nonrotating black rings in five 
dimensions, in particular the presence of a conical singularity. 

1 Introduction 



Perhaps the most interesting development of the last decade in black hole physics was Emparan and Reall's 
discovery of the black ring solution in cZ = 5 spacetime dimensions [1] [2] . This asymptotically flat solution 
of the Einstein equations has a horizon with topology S'^ x S^, while the Myers-Perry black hole |3j has a 

■ horizon topology S^. This made clear that a number of well known results in d = 4 gravity do not have a 
04 ■ simple extension to higher dimensions. 

■ The rapid progress following the discovery in (TJ [5] provided a rather extensive picture of the solutions' 
landscape for the five dimensional case [3], including configurations with abelian matter fields [5]-[TU]. Phys- 

■"sj" ' ically interesting solutions describing superposed black objects (black saturns [TT], bycicling black rings [12], 

, di-rings [T3j, concentric rings [14j) were also constructed. In all cases, the rotation provides the centrifugal 

repulsion that allows a regular solution to exist. In the static limit, all known d = 5 asymptotically flat 
solutions with a nonspherical horizon topology possess a conical singularity or other pathologies [2] , [5] , [10] . 

Despite the presence of several partial results in the literature, the d > 5 case remains largely unexplored 
(for a review, see [H]). Moreover, it becomes clear that as the dimension increases, the phase structure of 
^ • the solutions becomes increasingly intricate and diverse. The main obstacle stopping the progress in this 

5^ , field seems to be the absence of closed form solutions (apart from the Myers-Perry black holes), which were 

very useful in d = 5. Although the discovery in [TJ[2] was based on ^educated guesswork\ (and knowledge of 
the four dimensional C-metric), several systematic methods have been developed to construct d = 5 black 
ring solutions and their generalizations. Since no general framework seems to exist for d > 5, the issue of 
constructing black objects with a nonspherical horizon topology was considered by using approximations, in 
particular the method of asymptotic expansions [16j.[17j. The central assumption here is that some black 
objects can be approximated by a certain very thin black brane curved into a given shape. However, this 
method has limitations; for example it is supposed to work only if the length scales involved are widely 
separated. 

Perhaps a more straightforward approach to this problem would be to find such configurations numer- 
ically, as solutions of partial differential equations with suitable boundary conditions. This paper aims at 
a first step in this direction, since we propose a framework for a special class of solutions. As the simplest 
example of d > 5 black objects with a nonstandard topology of the event horizon, we present numerical 
evidence for the existence of static solutions with a S*^ x S'''^"' topology of the event horizon in d = 6 and 
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7 dimensions. These solutions share most of the features of the d = 5 static black ring, which is known in 
closed form 2 . In particular, they possess a conical singularity and approach the Schwarzschild-Tangherlini 
metric in a certain limit. 



2 A general ansatz and particular solutions 

2.1 The equations 

We consider the following metric ansatz: 

ds^ = -fo{r, z)dt^ + fi{r, z){dr^ + dz^) + ^r, z)d'^/ + h{r, z)dnl^^, (2.1) 

where dfl^_^ is the unit metric on S"^^^ and Q < r < oo, — oo < z < oo and ip is an angular coordinate, with 
< V" ^ 27r in the asymptotic region. Thus these coordinates have a rectangular boundary and are suitable 
for numerics. 

An appropriate combination of the Einstein equations, Gj =0, + = 0, = 0, and Gj] = (with 
GJ;^ the Einstein tensor), yields the following set of equations for the functions /i, /2, /a and /o 

VVo - ^(V/o)^ + :^(V/o) • (V/2) + ^^(V/o) • (V/3) = 0, 

VVi - j(^hf - iJ-m - 5)j||(V/3f - jI'^CV/o) ■ (VA) 

- %fi(v/.) . ,v/3) - 4^(vA)°(v/3, + (i^iMf = „. 

^/0/3 ^j2/3 J3 

VV2 - ^(V/2)' + ^(V/o) • (V/2) + ^^(V/2) • (V/3) = 0, (2.2) 
VV3 + ^^(V/3)^ + rV(V/o) • (V/3) + rV(V/2) • (V/3) - 2(d - 5)/i = 0, 

^/3 ^/O ^/2 

where we define (VU) ■ {VV) = drUdrV + d^Ud^V, V^U = dfU + dlU. 

The remaining Einstein equations G!! = 0, GJ^ — G^ = yield two constraints. Following [TSJ, we note 
that setting Gl ^ G^ ^ Gl + Gl ^ o'in V^G^'' = and V^G^^ = 0, we obtain the Cauchy-Riemann 
relations 



{V^gGl) + dr (^^g\{Gl - Gl)^ = 0, {V~gG:) - {y^g\{Gl ~ G^)) ^ 0. 



(2.3) 



Thus the weighted constraints satisfy Laplace equations, and the constraints are fulfilled, when one of them 
is satisfied on the boundary and the other at a single point [TH] . 



2.2 (/ = 5 

For d = h, this is the standard coordinate system used to study static, axisymmetric solutions of the vacuum 
Einstein equations, in which case the sphere f2d_4 reduces to a single angular coordinate (with < < 27r). 
The above equations present in this case a variety of physically interesting solutions which can be uniquely 
characterized by the boundary conditions on the z— axis, known as the rod-structure. In this approach, the 
z— axis is divided into N intervals (called rods of the solution), [—00, zi], [zi, Z2],. ■ ■ , [z^^i, 00]. A necessary 
condition for a regular solution is that only one of the functions /o(0, z), /2(0, z), /3(0, z) becomes zero for 
a given rod, except for isolated points between the intervals. For the static case discussed here, a horizon 
corresponds to a timelike rod where /o(0, z) = while lim^^o foif, z)/r'^ > 0. There are also spacelike rods 
corresponding to compact directions specified by the conditions fa{0,z) — 0, linir^o fa{r, z) /r'^ > 0, with 
a — 2,3. A semi-infinite spacelike rod corresponds to an axis of rotation, the associated coordinate being a 
rotation angle. 
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One of the main advantages of this approach is that the topology of the horizon is automaticaUy imposed 
by the rod structure. For example, the rod structure of a d = 5 static black ring solution consists of a 
semi-infinite space-like rod [—00, zi] in the -^-direction (thus /2(0, 2;) = there), a finite time-like rod [zi, Z2] 
(/o(0, z) = 0), a second (and finite) space-like rod [z2, Z3] in the ^'-direction, where /2(0, z) = again, and a 
semi-infinite space-like rod [23,00] (/3(0, z) = 0) in the (^-direction (see Figure fb). The metric functions 
of the static black rin^ are given by [5] , [T5] 

. _ i?2 + 6 . _ {Ri+^i+R2- 6)((1 - c)Ri + (1 + c)R2 + 2ci?3) 

8(1 + c)i?ii?2i?3 

^ (i?2-e2)(i?3+6) 



Ri-^1 



/a = i?3 - 6 , (2.4) 



where 



= R'l^Jr^+if and zi = -a, Z2 = a, Z3 = fe, (2.5) 



a and 6 being two positive constants, with c ~ a/b < 1. The mass, event horizon area and Hawking 
temperature of the d — 5 static black ring are: 



(with V3 — 27r^ the area of the three-sphere and G the Newton constant). Although this solution is asymp- 
totically flal|, it contains a conical excess angle S for the finite ■i/'-rod 



2.3 The Schwarzschild-Tangherlini black hole in d > 5 dimensions 

In principle, some d > 5 black objects can also be constructed by using a similar construction based on 
imposing a rod structure[f|. The starting point is the observation that the Minkowski spacetime is recovered 
within the framework (12.11) for 



/o = l, /i--^==, /2 = + z^ + z, h - v/r2 + z^ - z, (2.8) 
2vf + z 

(note that these expressiont|f| do not depend on the spacetime dimension d) . One can see that /2 vanishes 
for r = 0, z < 0, and /a for r = 0, z > 0, which, in the language of the Weyl formalism, corresponds to two 
semi-infinite rods [— oo,0] and [0, cxd]. 

Moreover, black holes with a spherical event horizon can also be written within the ansatz (j2.1|) as proven 
by the d > 5 Schwarzschild-Tangherlini solution. A straightforward but cumbersome computation leads to 
the following expression for the metric functions in this case 

/^(<i-3)/2 _ 1\ ^ (i;(d-3)/2 _^ l)4/(d-3) 

" i t,('i-3)/2 + 1 ) ' 2 (t-;^-!)^ I ,.2 (";^W T' ^^-^^ 

1 ('„('i-3)/2 , 1 A4/(d-3) 1 C„(d-3)/2 , i ^4/(d-3) 

/2 = ^ {c{v' + 1) + 2z.), h = ^ {c{v^ + 1)' 2zv), 



^Note that the function /i(0, z) behaves as l/|z — 2i| as z — > Zi. 

^The d = 5 static black ring solution in ^ admits an alternative interpretation as a ring sitting on the rim of a membrane 
that extends to infinity (this is found by requiring that the periodicity of tp is 2tt on the finite i/)-rod). However, the asymptotic 
metric is a deficit membrane in this case. 

^This possibility was already mentioned in Ref. 20 , which, to our knowledge, is the only attempt in the literature to 
numerically construct asymptotically Minkowski solutions with a nonspherical horizon topology. 

*The coordinate transformation r = sin 26, z = ^i?^ cos 29 (with < R < 00, < < ■7t/2), leads to a more common 
form of the flat spacetime metric: ds^ = —dt^ + dH? + B?{d9'^ + cos^ ddip'^ + s\v? 9dQ?^_^). 



3 



(a) 



(b) 
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t 



Q. Q. 



V \\i 

-c c -a a b 

Figure 1: The rod structure of the solutions is shown for the Schwarzschild-Tangherlini solution (left) and the black 
hole with 5*^ x 5"*~* topology of the horizon (right). The thin lines denote the z— axis and the thick lines denote the 
rods. Q stands for all directions on the sphere. 

where 

v^^(^r^ + z^ + V + V2y/{r^ + z'^f + c^ir'^ - z"^) + (r2 + z'^)V^ ^'"^ , 
and V = ^c^ + 2c2(r2 - z2) + (r2 + z2)2, 

with c an arbitrary positive constanl[f|. 

For any c? > 4, the rod structure of the Schwarzschild-Tangherhni black hole consists of a semi-infinite 
space-like rod [— oo, — c] (with 72(0,2)) = there), a finite time-like rod [— c, c] (/o(0, z) — 0) and a semi- 
infinite space-like rod [c, cx)] (with a vanishing /3(0, z)) in the f2-direction (see Figure la). The basic properties 
of the Schwarzschild-Tangherhni black hole can easily be rederived within this metric ansatz. Again, requiring 
the absence of a conical singularity imposes a periodicity 27r for the coordinate '!/'■ 

3 Black holes with x S'^~^ topology of the event horizon 

It is natural to conjecture that the ansatz (j2.ip can be used to construct d > 5 solutions with a nonspherical 
horizon topology. Perhaps the simplest case is provided by black holes with a topology of the horizon which 
is 5*2 X S'^^'^. The rod structure of these solutions is similar to the case of d = 5 static black rings, with the 
angle there replaced by the sphere, as shown in Figure lb. 

3.1 The rod structure and boundary conditions 

The new solutions in this paper have first a semi-infinite space-like rod [—00, —a] with the following expan- 
sio of the metric function as r — > (except for z — > ±a, 5): 

/o(r,z) = /ooW+rVo2W + ..., /i(r,2) = /io(z)+r2/i2W + ..., (3.1) 
/2(r, z) = r2/22(^) + rV24(2) + • ■ ■ , h{r. z) = h^{z) + h2{z) + . . . , 

the constraint equation = implying hm,-^o ''^/i//2 — ci- Next, there is a finite time-like rod [—a, a] 
corresponding to the horizon, with 

/o(r,z) =r2/o2(2)+rVo4W + ..., /i(r, z) = /io(z) + r2/i2 W + . . . , (3.2) 
/2(r, z) = /2o(z) + r'h2{z) + . . . , ^(r, z) = h^{z) + hi{z) + . . . , 

^ One can show that for d = 5 these expressions reduce to those in [2] . 

®In all relations below, the functions fik{z) arc solutions of a complicated set of nonlinear second order differential equations. 
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(with linir— >o ''^/i/Zo — const.), and a finite space-like rod [a, 6] along the ■(/'-direction, where the metric 
functions have an expression similar to p.ip : the constraint equation — implies there lim,.^o ^^/i//2 = 
C2, where ci, C2 are arbitrary positive constants. 

The case of the fi— rod {i.e. [b, oo]) is more involved. The expansion here is 

foir,z)=fooiz)+r\f22{z) + ..., /i(r, z) = /aa W + rVi2(^) + . . . , (3.3) 
/2(r, z) = /2o(z) + r^f22{z) + ..., fsir, z) = h2{z) + + . . . , 

Here the relation limr^o ''^/i/Zs = 1 follows, which is a feature of the d > 5 case. In five spacetime 
dimensions the value of this limit is not fixed. 

The obvious boundary conditions for large r, z is that fi approach the Minkowski background functions 
In practice, we have found it convenient to take 



/. = f^F,, (3.4) 

where are some background functions, given by the metric functions of the d — h static black ring solution 
(|2.4p . The advantage of this approach is that fi will automatically satisfy the desired rod structure. 

The equations satisfied by Fi can easily be derived from (|2.2p . As for the boundary conditions, the 
relations (|3.1[) - (|3.3p together with the expressions (|2.4p of the background functions /f imply 

drFi\r=o = 0, for — oo < z < 6, 

a.Fo|,=o = 0, a,Fi|,=o = 0, c',F2|,=o = 0, i^i|.=o = i^3|r=o for b<z<^, (3.5) 

and i^i = 1 as r ^ CX3 or z ±oo. 

The constraint equation = implies F^jFi = const, on the V-'-rods. Now, to be consistent with the 
assumption of asymptotic flatness, one finds const. = 1 for —oo < z < —a. The value of this ratio for the 
finite rod with a < z < b results from numerics. 

3.2 Global quantities 

The horizon metric is given by 

da^ = fiiO, z)dz^ + f2i0, z)diP^ + fsiO, z)dnl_^, (3.6) 

with —a < z < a. Since the orbits of shrink to zero at —a and a while the area of S"*^"^ does not vanish 
anywhere, the topology of the horizon is 5*^ x S'^~'^. 
The event horizon area is given by 

dzJhf2ft^ = 2AV'Fd_4^^ / dz{b- zy''-'y'JF,F2Ft\ (3.7) 
-a ^ Va + b J -a ^ 

where Vd-4 is the area of the unit sphere S'''"^ and Atp the periodicity of the angular coordinate ip on the 
horizon. 

The Hawking temperature can be computed from the surface gravity or by requiring regularity on the 
Euclidean section 



T=—xr-^J^, (3.8) 
47ra V 2 V -Fi 

where the constraint equation = guarantees that the ratio Fq/Fi is constant on the event horizon. At 
infinity, the Minkowski background is approached, with Aip — I-k there. The mass M of the solutions can 
be read from the asymptotic expression for /p 

^ (d-2)T/d_2(r2 + z2)('i-3)/4 + ■ ■ ■ ■ (^-^^ 
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Figure 2: The profiles of the metric functions Fi for a typical d — 6 solution with a — 0.5, b = 1.7. 



The solutions satisfy also the Smarr law 



M 



d-2 
40(^-3) 



TA 



H 



(3.10) 



All solutions we have found have F^/Fi 7^ 1 + 2a/{b — a) for the finite ^— rod. Thus, for a < z < 6, the 
{z,tp) part of the metric describes a surface that is topologically with a conical singularity at one of the 
poles. This conical singularity prevents this black object to collapse to form a spherical black hole horizon. 
The value of the conical excess of -0 for a < z < b is 



2n 1 



b + a F2 
b-a \ Fi 



We have found it convenient to introduce the quantity 

6/{27t) 



S = 



l-<5/(27r)' 

which has a finite range, — 1 < ^ < 0, and measures the 'relative angular excess'. 



(3.11) 



(3.12) 
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Figure 3: The metric functions Fi are shown as a function of z for several values of r. The relevant parameters here 
are d = 6, a — 1, b — 2.85. 



4 d = 6,7 numerical solutions 

We solve the resulting set of four coupled non-linear elliptic partial differential equations numerically, sub- 
ject to the above boundary conditions. All numerical calculations are performed by using the program 
FIDISOL/CADSOL [21j, which uses a Newton- Raphson method. 

First, one introduces the new compactified coordinates x ~ r/(l + r), u ^ arctan(z). The equations for 
Fi are then discretized on a non-equidistant grid in x and u. Typical grids used have sizes 80 x 160, covering 
the integration region < x < 1 and — 7r/2 < m < 7r/2. (See [21] and [22] for further details and examples 
for the numerical procedure.) In this scheme, the input parameters are the positions of the rods fixed by a 
and b and the value d of the spacetime dimension. To obtain solutions with S*^ x S''^"^ horizon topology, 
one starts with the c? = 5 solution as initial guess {i.e. Fi = 1) and increases the value of d slowly. The 
iterations converge, and, in principle, repeating the procedure one obtains in this way solutions for arbitrary 
d. The physical values of d are integers. We have studied solutions in d = 6, 7 dimensions in a systematic 
way. Solutions with d > 7 are also likely to exist; however, their study may require a different numerical 
method. 

The problem has two length scales a and 6, roughly corresponding to the event horizon radius and the 
radius of the round r2-sphere. It turns out that the above approach works well in the intermediate region 
of the ratio a/b, and fails to provide solutions with good accuracy if the two scales are widely separated or 
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Figure 4: (a): The scale free ratio Ah/M^"^ 2)/(d 3) shown as a function of the relative angular excess S as defined 
by the relation (|3TT2)) . The value of the parameter Cd there is Cd = (d - 2)/{l6TT)'^''~'^^^'-''~^^Vjlif~^^ and has been 
choosen such that the point (1,-1) on the plot corresponds to the Schwarzschild-Tangerlini black hole, (b): The 
relative angular excess 6 is shown as a function of the ratio between the two length scales a/b. 

very close to each other (i.e. a/b — > or a/b 1). FIDISOL/CADSOL provides also an error estimate for 
each unknown function. For d — 6, the typical numerical error for the functions is estimated to be lower 
than IQ-^. 

The Hawking temperature, event horizon area and conical excess are encoded in the values of the functions 
Fi, at r = 0. The mass parameter M is computed from the asymptotic form (j3.9p of the metric function 
gtt = —fo, the Smarr relation (|3.10p being used to verify the accuracy of the solutions. 

The functions Fi change smoothly with the two parameters a and b. Typical profiles of the solutions 
are presented in Figures 2, 3. One can see that the functions Fi are smooth outside of the z— axis and 
show no sign of a singular behaviour. The crucial point here is that the divergent behaviour of fi was 
already subtracted by the background functions fi- We have verified that the Kretschmann scalar stays 
finite everywhere, in particular at r = 0. 

By fixing the position of the horizon and varying the length of the finite -i/;— rod, we have generated a 
branch of d = 6, 7 black objects with S*^ x S'^^'^ horizon topology. The picture we have found is very similar 
to that valid for d — 5 static black rings. First, all solutions have a conical excess S on the finite ^— rod. 
Moreover, on the horizon, the area of the round S'''~'*-sphere is maximum for z = —a and minimum at z = a, 
which shows that Ai/j = 2tt there. 

In terms of the quantity S as defined in (|3.12p . these black objects smoothly interpolate between two limits 
(although these regions of the parameter space are difficult to approach numerically). First, as a/5 —> 1, one 
finds J — > — 1 {i.e. the conical excess 5 — > — oo) and the Schwarzschild-Tangherlini metric is approached (the 
finite -0— rod vanishes). As the second ip—rod extends to infinity {a/b 0), the radius on the horizon of the 
round S''^~''-sphere increases and asymptotically it becomes a {d — 4)-plane, while (5 — > 0. Here one expects 
to recover, after a suitable rescaling, the four dimensional Schwarzschild black hole uplifted to d dimension^ 
{i.e. a black (d — 4)-brane). These features are illustrated in Figure 4. For completeness, we have shown also 
the corresponding data for the d — 5 static black rings. 

^For the static d = 5 black ring solution, this limit is found by taking r v^26f, z —> •/2bz, ip —* Cp/y/Tb, followed by b ^ oo. 
This results in the Schwarzschild black string solution. 
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5 Further remarks 



In the absence of analytical methods to construct c? > 5 black hole solutions with nonspherical horizon 
topology, a numerical approach of this problem seems to be a reasonable task. In this work we have 
presented such a construction for d = 6,7 static black objects with S*^ x S'^^'^ topology of the horizon, as a 
first step before approaching more complex situations. The existence of the d — 6 solution is not a surprise, 
since the results in [23 show that this topology is one of the few allowed for black holes in six dimensions. 
However, given the presence of a conical excess angle, the solutions we have found are presumably unstable 
and their physical relevance is obscure. 

The technique proposed in this paper can easily be extended for other types of d > 5 static black objects 
{e.g. a superpositions of a Schwarzschild-Tangerlini black hole and a configurations with a S'^ x S'^^'^ topology 
of the horizon). Also, in principle, the inclusion of rotation in this scheme is straightforward. For example, 
similar to the d = 5 case, the conical singularity of the solutions with 5*^ x S''^"^ topology of the horizon 
could presumably be eliminated if the fi-sphere would rotate. However, this leads to a difficult numerical 
problem, since the equations would depend on at least three variables. Another possible direction would be 
to consider more complicated versions of the static ansatz (|2.1[) in d > 6 {e.g. replace the angular direction 
"0 with a sphere S*"). In the absence of rotation, one would expect such vacuum solutions to possess again 
some unphysical features. 

Moreover, the inclusion of some matter fields is unlikely to cure the conical singularity. This is the case 
for the Einstein-Maxwell-dilaton generalization of the solutions here, which extremize the action 

where F = dA. Starting with a vacumn configuration ()2.ip . a Harrison transformation (see e.g. [24]) leads 
to an EMd solution with line element 

ds2 ^ -/o(cosh2 /3 - sinh^ /3/o)-2"('*-3)di2 + (cosh^ (3 - sinh^ f3fo)^°'{fi{dr^ + dz^) + /ad^'^ + hd^ll_^), (5.2) 

and matter fields 

A^ = ^2{d-2)a \ 5^u = 0o - 2a(^ - 2)alog(cosh2 /3 - sinh^ /3/o), (5.3) 

cosh p — smh pjo 

where /3, (/)o are arbitrary real constants and a = {2a?'{d—2)-\-d—'i)~^ . Note that the Harrison transformation 
does not affect the rod structure of the metric. Thus, one finds in this way charged black objects with 
X S''"^ topology of the horizon. For d = 5, this is the solution presented in [5] by using a different 
coordinate system. 

Unfortunately, one can easily see that for both d — b and d = 6, 7, the line element (j5.2p has the same 
conical singularity as the vacuum seed solution. In principle, this singularity can be removed by "immersing" 
the solutions in a background gauge field, which requires to apply a second Harrison transformation (note 
that the resulting solutions would not be asymptotically flat; for d = 5, such a construction has been 
presented in [8]). The seed solution here is (|5.2[) . (|5.3p . with the electric field dualized, F = e^^"-'^ * F (with 
F = dB), and (j) = —4>. Then, the conical singularity vanishes for a critical value of the parameter in the 
second Harrison transformation. However, this last point requires knowledge of the explicit form of the 
magnetic potential B on the i/'-rods, which seems not possible for the d > 5 numerical solutions. 

We close this work by remarking that, on general grounds, a numerical approach works if the length scales 
involved are not widely separated, which is just the opposite of the approximate construction in [T6].|17j. 
Thus these methods are complementary. It would be interesting to construct solutions with S'^~^ x 
topology of the event horizon, which were already considered in [16j. Supposing that the static limit of these 
black rings can be constructed by using the ansatz (|2.ip , their rod structure can be read from Figure lb by 
interchanging the ip and 17-rods there {e.g. there is only one ip—iod for [&, oo], etc.). The boundary condi- 
tions in this case are still given by p.ip - (|3.2p . An unexpected feature here is that these configurations would 
possess no conical singularities. This difference to the static d = 5 black rings originates in the presence of 
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{d — 5) factors in the field equations l|2.2p ). This counter-intuitive result sheds doubt on the possibility to 
construct d > 5 static black rings with a regular horizon within the ansatz (|2.ip with a dependence on only 
two coordinates. 
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